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retains  recognition  specificity  through  rare,  molecule -unique  coincidences  between  IR  molecular  absorption  and  a 
line-tunable  C02  excitation  laser.  The  laser-induced  double  resonance  is  observed  as  a  modulated  THz  spectrum 
monitored  by  a  THz  transceiver.  As  an  example,  our  analysis  indicates  that  a  1  ppm  cloud  of  CH3F  100  m  thick  can 
be  detected  at  distances  up  to  1  km  using  this  technique. 

Keywords:  Infrared,  Terahertz,  Spectroscopy,  Remote  Sensing, 


1.  INTRODUCTION 

Double  resonance  (DR)  spectroscopy  is  a  form  of  wavelength  non-degenerate  pump-probe  spectroscopic  analysis. [1] 
For  molecular  DR  spectroscopy,  pulsed  infrared  (IR)  radiation  excites  gas  phase  molecules,  creating  a  non- 
equilibrium  distribution  of  population  in  quantized  rotational  (and  vibrational)  states  that  may  be  probed  by  terahertz 
(THz)  radiation.  Because  the  spacing  between  rotational  energy  levels  (-1  cm"1)  is  much  less  than  kT  (200  cm'1  at 
300K)  while  the  energies  of  excited  vibrational  levels  are  typically  several  kT  above  the  ground  state,  only  the 
rotational  manifold  in  the  ground  vibrational  state  will  be  heavily  populated.  When  IR  pump  radiation  is  tuned  to 
coincide  with  a  specific  ro- vibrational  molecular  transition,  the  pump  excites  molecules  from  one  of  the  heavily 
populated  rotational  states  in  the  ground  vibrational  level  to  one  of  the  sparsely  populated  rotational  states  in  an 
excited  vibrational  level.  The  dynamical  evolution  of  the  non-equilibrium  population  distributions  in  both  rotational 
manifolds  can  be  probed  with  terahertz  radiation  tuned  to  coincide  with  one  of  the  pure  rotational  transitions  in  the 
respective  manifold.  The  strongest  DR  signals  will  occur  for  those  rotational  transitions  containing  one  level 
directly  affected  by  the  IR  pump,  and  it  is  these  molecule -unique  signals  that  present  an  opportunity  for  trace  gas 
remote  sensing. 

In  order  to  examine  the  potential  of  this  technique  for  remote  sensing  at  atmospheric  pressures,  we  will 
quantitatively  compare  the  expected  strength  of  these  DR  signals  with  various  natural  and  system  noise  sources  to 
estimate  the  detectable  signal  to  noise  ratio  (SNR)  for  various  operational  scenarios.  This  manuscript  will  examine 
the  gas  methyl  fluoride  (CH3F)  because  it  has  been  extensively  characterized  at  low  pressures  in  order  to 
demonstrate  how  DR  spectroscopy  can  overcome  the  limitations  imposed  by  atmospheric  pressure  broadening. [2,3,4] 
An  introduction  to  low  pressure  spectroscopic  techniques  will  be  followed  by  estimates  of  THz  DR  signatures  and 
their  detectability  in  a  representative  operational  scenario. 

2.  SINGLE  RESONANCE  SPECTROSCOPY 

Low-pressure  (<100  mTorr)  rotational  spectroscopy  provides  an  exquisitely  precise  technique  for  molecular 
recognition.  The  spectra,  which  can  be  routinely  measured  with  an  accuracy  of  1 : 1 07,  depend  sensitively  on  the 
shape  (i.e.  moments  of  inertia  and  atomic  masses)  of  the  molecule. [5,6]  For  example,  a  typical  low-pressure 
rotational  spectrum  of  the  prolate  symmetric  top  molecule  methyl  fluoride  (12CH3F)  is  depicted  in  Figure  1. 
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Figure  1:  (Left)  12CH3F  molecule.  (Right)  12CH3F  ground  vibrational  state,  J  =  12— ►  13,  K  =  0,1,  2,...  ,12  rotational  spectrum  is 

shown. 

For  comparison,  the  spectrum  of  the  isotopic  isomer  13CH3F  (1%  isotopic  abundance)  occurs  at  645  GHz,  well 
separated  from  the  spectrum  of  12CH3F.[7]  The  lines  in  Figure  1  are  resolvable  because  they  were  obtained  at  low 
pressure  (10  mTorr)  where  linewidths  are  <  1  MHz;  however,  at  atmospheric  pressure  (760  Torr)  specificity  is  lost 
because  atmospheric  pressure  broadening  produces  linewidths  of  ~  5  GHz.[8] 

Although  not  as  precise  as  THz  rotational  spectroscopy,  IR  ro-vibrational  spectroscopy  is  widely  used  for  molecular 
recognition  because  of  the  availability  of  IR  spectrometers  and  the  transparency  of  the  atmosphere  in  the  well- 
known  IR  windows. [9]  The  IR  spectrum,  which  can  be  routinely  measured  with  an  accuracy  of  1 : 1 05,  depends  on 
the  type  of  bond  that  is  vibrating. [9]  For  example,  the  infrared  spectrum  of  the  C-F  stretch  mode  in  12CH3F  is 
depicted  in  Figure  2. 
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Figure  2:  Fourier  transform  IR  spectra  of  C-F  stretch  (V3  =  1  level)  in  12CH3F. [10,11] 


Ro-vibrational  transitions  occur  in  three  distinct  types,  determined  by  the  whether  the  rotational  state  changed  by  AJ 
=  -1,0,  or  l.[9]  These  are  called  P-branch,  Q-branch,  and  R-branch  transitions,  respectively,  and  are  distinguishable 
as  three  separate  bands  from  right  to  left  in  Figure  2.  IR  transitions  are  labeled  by  the  ground  rotational  state 
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involved  and  the  type  of  IR  transition,  e.g.  P(J)  for  C02  or  Q(J,K)  for  symmetric  top  molecules.  Note  that  Q-branch 
transitions  are  not  as  well  resolved  as  the  P-  and  R-branch  transitions. 


3.  DOUBLE  RESONANCE  SPECTROSCOPY 

Unlike  traditional  single  resonance  spectroscopy,  which  measures  the  absorption  of  a  transition  based  on  the 
equilibrium  population  difference  between  the  upper  and  lower  states,  DR  spectroscopy  measures  the  change  in  the 
population  difference  induced  by  the  pump.  The  variation  in  absorption  of  a  rotational  transition  can  be  monitored 
by  measuring  the  change  in  coincident  terahertz  radiation  transmitted  through  the  molecular  sample.  If  the  pump 
excitation  is  strong  enough,  a  rotational  population  inversion  can  be  induced.  This  is  the  basis  for  low  pressure 
optically  pumped  far-infrared  (OPFIR)  lasers. [12] 

OPFIR  lasers  illustrate  the  importance  of  a  frequency  coincidence  between  the  pump  (a  C02  laser)  and  a  molecular 
ro-vibrational  transition.  Although  a  typical  molecule  will  have  thousands  of  ro -vibrational  transitions,  it  is  rare  if 
any  of  those  transitions  coincide  with  any  line  from  a  C02  laser.  For  these  coincidences  to  occur,  two  criteria  must 
be  met.  First,  there  must  be  a  vibrational  level  within  the  tuning  range  of  a  C02  laser:  between  850cm"1  and  1150 
cm'1.  This  criterion  becomes  easier  to  meet  the  heavier  and  more  complex  the  molecule  becomes. [13]  Second,  a 
specific  molecular  ro-vibrational  transition  must  coincide  with  a  specific  C02  laser  line  to  within  the  tuning  range  of 
the  laser  (~  100  MHz)  and  the  linewidth  of  the  transition  (~  30  MHz  at  low  pressure). 

When  the  IR  spectrum  for  the  C-F  stretch  in  12CH3F  (red)  overlays  the  available  C02  laser  lines  (green),  the 
possibility  of  coincidences  between  the  two  emerges  (Figure  3).  Although  it  appears  that  there  should  be  many 
naturally  occurring  coincidences,  there  are  only  two  close  ones  at  low  pressures:  the  9P(20)  C02  laser  line  can 
excite  the  Q(12,l)  or  Q(12,2)  transitions.  [14]  Likewise,  only  one  close  coincidence  occurs  at  low  pressure  for 
13CH3F:  the  9P(32)  C02  laser  line  excites  the  R(4,3)  transition. 
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Figure  3:  C02  laser  lines  are  depicted  in  green,  spectra  of  12CH3F  is  depicted  in  red. 


4.  COINCIDENCE  MODELING  and  DOUBLE  RESONANCE  SIGNATURES 

One  of  the  tremendous  advantages  of  performing  DR  spectroscopy  at  atmospheric  pressure  is  that  the  ro-vibrational 
transitions  of  the  trace  gas  (e.g.  CH3F)  broaden  from  ~30  MHz  to  ~  5  GHz,  due  to  rapid  collisions  with  atmospheric 
N2  and  02,  dramatically  relaxing  the  coincidence  requirement.  To  estimate  how  many  more  coincidences  occur,  the 
ro-vibrational  transition  frequencies  are  calculated  for  all  thermally  populated  rotational  states.  The  frequencies  of 


3 


these  transitions  are  then  compared  to  the  catalog  of  available  C02  laser  frequencies  in  order  to  identify 
coincidences  that  are  within  5  GHz.  In  this  way,  many  new  coincidences  are  found,  as  illustrated  in  a  portion  of  the 
12CH3F  Q-branch  spectra  is  depicted  in  Figure  4. 
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Figure  4:  Predicted  coincidences  (green)  between  a  trace  amount  of  12CH3F  [green,  Q(J,K)]  in  the  atmosphere  and  C02  laser 

lines  [red,  9P(J)]. 

When  a  ro-vibrational  transition  is  excited  by  a  coincident  C02  laser  line,  four  rotational  transitions  are  directly 
connected  to  the  two  ro-vibrational  levels  (Figure  5).  Regardless  of  whether  the  rotational  transition  is  in  the  excited 
or  ground  vibrational  level,  the  change  in  absorption  of  each  of  the  four  transitions  will  be  approximately  the  same: 
the  pump  removes  the  same  number  of  molecules  from  the  “ground  state”  as  it  adds  to  the  “excited  state”.  Of 
course,  the  sign  of  the  change  will  be  different,  depending  on  whether  the  transition  is  above  or  below  the  pumped 
ro-vibrational  level.  Two  of  the  rotational  transitions  will  exhibit  a  pump-induced  increase  in  absorption,  and  two 
will  exhibit  a  decrease  in  absorption  or  perhaps  even  emission.  The  separation  between  the  components  is  given  by 
the  spacing  of  the  rotational  levels,  which  is  roughly  2B,  where  the  rotational  constant  B  ~  25  GHz  for  CH3F.[  6,15] 
Note  that  this  50  GHz  of  separation  allows  the  DR  signature  to  be  well  resolved,  even  at  atmospheric  pressure  where 
the  linewidths  are  ~5  GHz. 

From  Figure  5,  it  is  clear  that  the  resulting  THz  DR  signature  depends  on  the  type  of  ro-vibrational  transition 
involved:  P-,  Q-,  or  R-branch.  P-  and  R-branch  transitions  induce  opposite  DR  spectra  whose  two  middle  features 
are  of  the  same  sign  and  nearly  overlap.  Q-branch  transitions,  on  the  other  hand,  create  two  doublets  of  opposite 
sign  whose  components  nearly  overlap.  These  near-overlaps  arise  from  the  difference  between  the  rotational 
constants  in  the  ground  and  excited  vibrational  levels,  differences  that  shrink  with  increasing  molecular  mass.[6] 
These  DR  THz  spectra  are  modulated  signals  that  are  only  detectable  when  the  pulsed  laser  is  exciting  the  gas  and 
before  collisions  redistribute  the  pumped  molecules  throughout  the  rotational  manifold.  The  timescale  for  collisional 
redistribution  is  determined  by  the  mean  collision  time  at  atmospheric  pressure,  namely  t  =  1/(2*7t*2.5  GHz)  ~  100 
ps.[8]  Therefore,  the  DR  signal  that  uniquely  identifies  the  molecule  must  be  excited  and  measured  during  this  time. 
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Figure  5.  (Left)  Notional  energy  diagram  of  the  ro -vibrational  states  affected  by  the  pump  laser  for  the  three  IR  transition  types 
(P-,  Q-,  and  R-branch),  respectively.  (Center)  A  graphical  depiction  of  the  corresponding  pump -induced  DR  signal.  (Right) 

Estimated  frequency  dependent  DR  signal  for  each  case. 

5.  SPECIFICITY 


The  attributes  of  rare,  molecule -specific  pump  coincidences  and  well-resolved  DR  THz  spectra  make  DR 
spectroscopy  a  viable  technique  for  detection  of  trace  gases  in  the  atmosphere  in  ways  exceedingly  challenging  for 
more  traditional  single  resonance  THz  spectroscopy.  A  two-dimensional  specificity  matrix  can  be  compiled  that 
indicates  the  specific  combination  of  IR  and  THz  frequencies  that  uniquely  identify  a  molecule,  even  distinguishing 
isotopic  isomers.  For  example,  12CH3F  is  excited  by  the  9P(20)  line  to  produce  Q-branch  DR  THz  signatures  with 
features  at  604.3,  612.4,  654.6  and  663.4  GHz.[7]  13CH3F  is  excited  by  the  9P(32)  line  to  produce  R-branch  DR  THz 
signatures  with  features  at  198.8,  248.6,  245.4,  and  294.4  GHz.[16]  This  is  an  enormous  difference  for  two  molecules 
that  have  identical  structures  and  only  differ  in  mass  by  1  amu. 
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Figure  6:  The  2D  specificity  matrix  for  12CH3F  double  resonance  remote  sensing.  The  most  practical  DR  spectra  for  remote 

sensing  are  highlighted  in  green,  the  least  in  red. 


As  the  pressure  increases,  the  number  of  coincidences  increases,  so  that  the  most  convenient  combination  of  pump 
(IR)  and  probe  (THz)  frequencies  may  be  selected.  Figure  6  graphically  illustrates  the  specificity  matrix  of  DR 
signatures  for  each  combination  of  laser  line  and  THz  frequency  for  12CH3F.  New  coincidences  that  are  predicted  to 
occur  at  atmospheric  pressure  are  presented  above  the  well-known  low  pressure  coincidence  induced  by  the  9P(20) 
line.  The  strength  of  the  DR  THz  signatures  depends  strongly  on  IR  absorption  strengths,  laser  line  power, 
atmospheric  absorption,  and  THz  source  availability.  The  potential  utility  of  each  coincidence  is  ascertained  next. 


6.  DOUBLE  RESONANCE  ABSORPTION  STRENGTH 


To  estimate  the  strength  of  the  pump-induced  DR  THz  signature,  consider  first  the  single  resonance,  frequency- 
integrated  THz  absorption  coefficient  ay  for  any  of  the  four  transitions  involved: 


a, 
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=  D - lL 


3hc 


(i) 


where  i  and  j  are  adjacent  rotational  states  in  the  same  vibrational  level,  V y  and  (nrnj)  are  the  (THz)  rotational 
transition  frequency  and  equilibrium  population  difference  between  them,  and  ju'y  is  the  dipole  matrix  element  for 
that  rotational  transition. [8]  In  order  to  calculate  the  peak  absorption  at  a  given  frequency,  the  integrated  absorption 
coefficient  must  be  multiplied  by  an  area-normalized  Lorentzian  lineshape  function,  S(v,v0),  since  the  transition  is 
pressure  broadened  at  atmospheric  pressure. 


Under  thermal  equilibrium  conditions,  ^  and  rij  are  determined  though  Boltzmann  distributions  and  are  dependent 
upon  the  total  number  of  molecules  present  in  the  molecular  gas  in  a  given  volume,  n0  =  P/RT ,  where  P  is  the  partial 
pressure  and  T  is  the  temperature  of  the  trace  gas.  For  a  trace  amount  of  target  molecular  gas  in  an  atmosphere  the 
peak  absorption  is  reduced  by  pressure  broadening  by  a  dilution  factor,  D,  which  is  simply  the  ratio  of  the  partial 
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pressure  of  the  gas  to  atmospheric  pressure.  For  example,  the  peak  absorption  for  1  ppm  of  CH3F  is  reduced  by  D  = 
10"6  over  the  peak  absorption  of  pure  CH3F  of  the  same  (partial)  pressure. 

Turning  now  to  calculate  the  DR-induced  change  in  absorption,  assume  that  the  pump  pulse  is  comparable  to  the 
atmospheric  collision  time  so  that  there  is  no  collisional  thermalization  of  excited  molecules.  The  pump  rate  R^  of  a 
particular  ro-vibrational  transition  i  — jy  by  the  C02  laser  with  power  p  per  pulse  and  beam  radius  r  is  given  by 


(2) 


where  V 0  and  //L  are  the  frequency  and  ro-vibrational  matrix  element  connecting  the  ground  and  excited  vibrational 
state,  and  the  pressure -broadened  (Lorentzian),  area-normalized  line  shape  function  S{vp,  \?0)  is  included  to  account 
for  the  reduced  pumping  efficiency  when  vp  ±  i/0.[17,18,8] 


From  equation  2  it  is  evident  that  as  the  C02  pump  irradiance  {p/nr2)  increases,  the  number  of  molecules  that  can  be 
moved  per  second  from  the  ground  vibrational  state  to  the  excited  vibrational  state  increases.  The  number  of 
molecules  excited  by  the  pump  np  is  then  simply  the  product  of  R ij9  the  laser  pulse  width  x,  and  the  number  of 
molecules  in  ground  rotational  state  nh 


np  =  t  tit  =  Ry  zfi  nQ,  (3) 

where  f  is  the  fraction  of  the  total  molecular  population  in  state  i.  Since  it  is  assumed  that  there  is  no  collisional 
redistribution  of  population  during  the  brief  laser  pulse,  the  modified  single  resonance  THz  absorption  may  be 
estimated  through  the  change  in  population  of  the  affected  states  using  (1): 


<j=D~ 


3hc 


Mu  («,-«,+ 'O, 


(4) 


where  the  sign  preceding  np  is  determined  by  whether  the  pump  increases  (+,  green  arrows  in  Fig.  5)  or  decreases  (-, 
yellow  arrows  in  Fig.  5)  the  absorption  of  the  rotational  transition  in  question. 

Consider  the  specific  case  of  12CH3F  in  which  the  C02  laser  excites  the  Q(12,2)  ro-vibrational  transition.  For  the  J  = 
12— >13  transition  in  the  ground  vibrational  level,  the  pump  removes  np  molecules  from  the  equilibrium  population 
n 72,  so  the  resulting  a’12-i3  is  reduced  (Fig.  7(a))  or  even  inverted  if  np  >  n12  -  n13  (Fig.  7(b)).  Likewise,  in  the  excited 
vibrational  level  np  molecules  are  added  to  the  J  =  12  rotational  state,  causing  an  increased  THz  absorption  a’  12_13,  as 
shown  in  Figure  7(c). 
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Figure  7:  Equilibrium  and  pumped  THz  rotational  spectra  for  representative  transitions  in  12CH3F  at  a  lppm  concentration,  for 
which  the  Q(12,2)  transition  is  excited  by  a  9P(20)  C02  laser  line.  (Left)  The  J  =  12  — >13,  K  =  2  transition  in  the  ground 
vibrational  level  for  100  mJ  of  C02  laser  power  in  r  =  100  cm  beam  and  a  100  ps  pulse.  (Center)  The  same  transition  with  1.0  J  of 
C02  laser  power.  (Right)  The  same  transition  and  pump  parameters  as  in  (Center)  but  for  the  J  =  12 — >13,  K  =  2  transition  in  the 

excited  vibrational  level. 
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The  DR  THz  absorption  coefficient,  A  a ,  is  then  simply  the  difference  between  the  equilibrium  and  pumped 
absorption  coefficients: 


_  8 nzvru  | 

Aaij=a\-aij=+D-^—\M;. 


(5) 


This  A  ay  is  what  would  be  measured  when  remotely  sensing  trace  gases.  Notice  that  A  ay  depends  linearly  on  the 
pump  intensity  p/nr2.  More  subtly,  it  depends  on  the  vibrational  dipole  matrix  element  juvy,  which  generally  favors 
P-  and  R-  type  ro-vibrational  transitions  except  at  low  J.  It  is  important  to  recognize  that  during  the  short  duration  of 
the  pump  pulse  when  thermalizing  collisions  have  not  yet  occurred,  only  the  four  rotational  transitions  directly 
connected  to  the  pumped  ro-vibrational  transition  exhibit  any  THz  DR  signature.  In  fact,  they  have  the  same  A  a  if. 
Aa12-i3tGs  ~  -Aau_i2,GS  ~  AanG2.ES  ~  - Aa12.i3,ES •  For  all  other  rotational  transitions  A  ay  =  0.  When  these  four 
components  are  plotted  with  the  appropriate  lineshape  function  centered  on  their  respective  frequencies,  the  total  DR 
signature  for  this  molecule  emerges  (Figure  8). 
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Figure  8:  Constituent  THz  DR  spectra  in  the  excited  vibrational  level  (Top)  and  ground  vibrational  level  (Center),  as  well  as  the 

sum  DR  signature  that  would  be  observed  (Bottom). 
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One  final  point  must  be  made.  Because  of  the  broad  lineshapes  at  atmospheric  pressure,  the  pump  laser  actually 
coincides  with  most  if  not  all  2J+1  of  the  K  components,  making  each  of  the  four  features  in  the  DR  signature  the 
sum  of  2J+1  components.  For  example,  in  12CH3F,  the  9P(20)  laser  line  is  actually  coincident  with  ten  of  the 
Q(12,K)  ro-vibrational  transitions  (K  =  0-9)  at  atmospheric  pressure.  The  DR  signatures  for  each  of  the  constituent 
K  components,  and  the  summed  DR  signature,  are  shown  in  Figure  9. 


THz  Frequency  (GHz) 


Figure  9:  (Left)  The  individual  DR  signals  from  each  of  the  ten  K  components  of  excited  Q(12,K)  transition  contributing  to  the 
summed  12CH3F  DR  signature  (Right).  Note  the  change  in  vertical  scale  on  the  right. 


7.  DOUBLE  RESONANCE  REMOTE  SENSING 

Now  that  the  procedure  to  estimate  THz  DR  signatures  has  been  described,  we  may  use  this  calculation  to  illustrate 
how  double  resonance  spectroscopy  may  be  used  for  remote  sensing  of  trace  gases [19].  Suppose  the  objective  was  to 
detect  a  100  m  thick  cloud  of  1  ppm  methyl  fluoride  at  a  distance  of  1  km.  Assume  the  IR  pump  beam  and  THz 
probe  beam  will  be  configured  so  that  they  completely  overlap,  travel  in  a  horizontal  path  through  cloud,  then  reflect 
back  to  the  THz  transceiver  from  a  naturally  occurring  or  cooperative  background  reflector.  For  long-term 
monitoring,  it  is  acceptable  to  use  a  cooperative  retroreflective  material,  such  as  copper  with  a  reflectivity  of  99%, 
oriented  normal  to  the  transceiver  with  a  size  that  is  at  least  as  large  as  the  diffracted  beam  at  1km. 

Atmospheric  attenuation  of  THz  radiation  presents  a  severe  restriction  on  the  applicability  of  this  technique.  THz 
DR  signatures  are  only  measurable  within  one  of  the  four  atmospheric  “windows”  identified  in  Figure  10,  whose 
transparency  depends  sensitively  on  the  water  vapor  content  of  the  atmosphere.  The  ambient  atmospheric 
attenuation,  as  well  as  the  single  resonance  THz  absorption  by  the  trace  gas  cloud  (equation  1),  must  be  calculated. 
For  the  former,  we  used  the  clean  air  Millimeter-wave  Propagation  Model  (MPM)  developed  by  Liebe.[20] 
Contributions  from  molecular  oxygen  and  water  are  included  in  the  MPM  calculations  along  with  the  non-re  sonant 
dry  air  spectrum  and  the  water  vapor  continuum.  For  the  12CH3F  THz  DR  signature  occurring  at  663  GHz,  double 
pass  atmospheric  attenuation  is  20  dB/km  for  a  25 °F  dew  point,  barely  measurable,  while  it  is  an  unmeasurable  162 
dB/km  for  a  75  °F  dew  point.  Luckily,  the  relaxed  coincidence  requirement  made  possible  by  atmospheric  pressure 
broadening  increases  the  probability  that  a  given  molecule  will  have  a  THz  DR  signature  in  one  of  the  lower 
frequency,  higher  transparency  windows. 
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Figure  10:  Calculations  of  atmospheric  attenuation  of  THz  propagation  for  selected  dew  points.  Four  atmospheric  windows 

amenable  for  DR  spectroscopy  are  indicated. 

By  combining  the  DR  signature  calculations  and  the  MPM-based  THz  propagation  models,  the  scenario-dependent 
DR  signal  strength  received  at  the  THz  transceiver  may  be  estimated  as  a  function  of  IR  pump  power  and  THz 
transmitter  power.  For  these  calculations  assume  that  a  mode-locked  C02  TEA  laser  provides  ten  100  ps,  1  J 
micropulses  per  macropulse  with  a  10  Hz  macropulse  repetition  rate  and  that  a  continuous  wave  THz  source  emits 
100  pW.  The  frequency  dependent  DR  spectra  generated  based  on  these  parameters  for  12CH3F  is  shown  in  Figure  9. 

To  ascertain  whether  or  not  this  DR  signal  is  detectable,  an  estimate  of  the  system  noise  must  be  made.  Two 
otherwise  insidious  noise  sources  that  plague  other  remote  sensing  techniques  -  blackbody  radiation  and 
atmospheric  fluctuations  -  are  negligible  for  this  technique,  the  former  because  there  is  virtually  no  blackbody 
radiation  in  the  narrow  bandpass  of  the  THz  heterodyne  receiver,  the  latter  because  atmospheric  fluctuations  are 
frozen  on  the  100  psec  timescale.  Of  the  remaining  possible  sources  of  noise,  two  forms  of  detector  noise  dominate: 
thermal  noise  and  Townes  noise. [6]  Because  relatively  large  THz  powers  will  be  required,  Townes  noise  will 
dominate  all  others: 


PN  =A/kT(bB)05Pc  (6) 

where  T  is  the  noise  temperature  of  the  THz  heterodyne  receiver,  b  is  the  IF  bandwidth,  and  B  is  the  post  detection 
bandwidth.  For  a  THz  source  power  of  lOOpW  the  received  carrier  power,  after  propagation  and  retroreflection,  is 
Pc  =  4.4e-7  W.  For  a  receiver  with  T=  1273  K,  b  =  10  GHz  and  i?  =  10  GHz,  the  Townes  noise  is  ~2.8  xlO"8  W. 

Combining  these  estimates  of  the  12CH3F  THz  DR  signal  and  noise  for  a  0  °F  dew  point,  the  estimated  Aay  produces 
a  laser-induced  change  in  the  received  THz  power  of  ~1.3xl0"7  W  per  laser  pulse,  and  for  a  one  second  integration 
time  (100  micropulses)  the  estimated  SNR  is  an  easily  detectable  46.4.  This  analysis  shows  that  12CH3F  becomes 
undetectable  at  higher  dew  points  because  the  atmospheric  window  at  650  GHz  closes.  However,  the  isotopic 
isomer  13CH3F,  whose  double  resonance  signature  occurs  near  250  GHz,  remains  detectable  in  spite  of  its  rare,  1% 
natural  abundance  even  up  to  the  highest  relative  humidity  (Figure  11).[19]  This  analysis  confirms  that  it  is  much 
more  favorable  to  operate  in  the  lower  frequency  THz  windows  because  of  the  reduced  sensitivity  to  water  vapor 
absorption.  Our  analysis  also  finds  that  as  molecules  become  heavier  and  the  B  rotational  constant  decreases,  more 
of  the  DR  THz  signatures  occur  at  these  lower  frequencies.  Therefore,  if  a  light  molecule  like  CH3F  can  be  detected 
by  this  technique,  heavier  molecules  are  even  more  likely  to  be  detected. 
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Figure  11:  Estimates  of  the  detectability  of  THz  DR  signatures  as  a  function  of  atmospheric  water  vapor  content  for  the  natural 
abundances  of  both  isotopes  of  CH3F,  as  described  in  section  7.  Calculations  are  performed  at  the  optimal  THz  frequency  for 
both  isotopes:  12CH3F  at  663  GHz  and  13CH3F  at  246  GHz.  The  dashed  line  indicates  an  SNR  of  1. 


8.  CONCLUSIONS 

The  feasibility  of  IR/THz  DR  spectroscopy  for  remote  sensing  has  been  demonstrated  through  an  examination  of  a 
prototypical  molecule,  CH3F  in  a  1  ppm,  100  m  thick  cloud  at  a  range  of  1  km,  detected  in  plausible  scenarios  using 
available  or  constructible  hardware.  The  sensitivity  to  pump  intensity,  ro -vibrational  transition  type,  probe  power, 
atmospheric  water  vapor  content,  and  location  of  the  THz  DR  signature  within  an  atmospheric  window  has  been 
investigated  to  outline  the  general  applicability  of  this  technique.  These  insights  have  revealed  certain  advantages 
for  working  at  atmospheric  pressure  and  have  indicated  a  favorable  scaling  to  increasingly  heavy  and  complex 
molecules. 
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